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The question of the variation of (1) the direct and (2) the diffused 
radiation is next discussed. On comparing the curve representing the 
chemical intensity of diffused light with the curve of solar obscuration, it 
is found that the rate of diminution in chemical action exerted by the dif- 
fused light is up to a certain point greater than corresponds to the portion 
of eclipsed sun, whilst from this point up to totality the rate of diminution 
becomes less than corresponds to the progress of the eclipse. The same 
rapid diminution in the chemical action of the diffused daylight during the 
early periods of the eclipse was also observed at Jamkandi ; it is doubtless 
due to the dark body of the moon cutting off the light from the brightly 
illuminated portion of sky lying round the solar disk. 

The results of the observations at Catania are then compared with those 
made at Moita, near Lisbon, and communicated to the Society in 1870. 
This comparison shows a striking coincidence between the two sets of ob- 
servations. In each case it is seen that the relation between solar altitude 
and total chemical intensity is represented by a straight line, although the 
Catania observations slightly exceed, by a constant difference, those made 
at Moita in conformity with the slight difference in latitude, and with the 
fact that the former determinations were made at a greater elevation above 
the sea-level. 

The Catania observations further confirm the fact which we formerly 
announced, that for altitudes below 50° the amount of chemical action 
effected in the plane of the horizon by diffused daylight is greater than that 
exerted by direct radiation, and also that at altitudes below 1 0° direct sun- 
light is almost completely robbed of its chemically active rays. 

VII. "On the Calculation of Euler 5 s Constant." By J. W. L. 

Glaisher, B.A., F.B. A.S. Communicated by James Glaisher, 
F.B.S. Beceived June 6, 1871. 

The main object of the present communication is to correct some inac- 
curacies both of reasoning and calculation contained in two papers by Mr. 
Shanks, viz. " On the Extension of the Value of the Base of Napier's 
Logarithms ; of the Napierian Logarithms of 2, 3, 5, and J ; and of the 
Modulus of Briggs on the common System of Logarithms ; all to 205 
places of Decimals," in the Proc. Boy. Soc. vol. vi. p. 397; and "On 
the Calculation of the Numerical Value of Euler ? s Constant," in the Proc. 
Boy. Soc. vol. xv. p. 429 (1867). 

For the calculation of the constant Mr. Shanks has used (as, indeed, 
has every calculator who has computed the value of the constant during 
the present century) the semiconvergent series 

2 3 x 2qg 2cc 4 or 6 a? 

y being the constant, and B l3 B s , B B . . . Bernoulli's numbers. 
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Mr. Shanks obtains the value of y from this formula by making ^=10, 
20, 50, 100, 200, 500, and 1000, and remarks, as a curious coincidence, 
that the number of decimal places obtained from x being made equal to 
10, 20, 50, and 100 is nearly proportional to |/ 10, 8,20, »/50, and £/100. 
On p. 432 he gives a Table of the number of decimal places obtainable 
from the formula when x has the values 2, 5, 10, 20, 50, 100, 200, 500, 
and 1000, and from it draws the inference that "we may fairly infer that 
when n is increased in a geometrical ratio, the corresponding number of 
decimals obtained in the value of E increases only in something like an 
arithmetical one, and that probably from 50,000 to 100,000 terms in the 
Harmonic Progression would require to be summed in order to obtain 100 
places of decimals in the value of E, Euler's constant/' 

Algebraically, of course, y is independent of x in the formula (i), but 
arithmetically, since the series ultimately becomes divergent, the value of y 
is so far dependent on x that for a given value of x the series will only 
afford a certain number of decimal places. The number of decimal places 
directly obtainable is equal to the number of ciphers which precede the first 

significant figure in the value of the numerically least term of the series. 

1 
The wth term (considering only the terms after ^~ in (i), so that the 

nth term is the term involving B„) 

B 

2nx %l% 
which* since 

B J (1.2»3„.2n) / 1 1 \ 

is very hearly equal to 

2(1.2.3. . . (2^-1)} . 

(2itx) 2 « ; 

so that the ratio of the nth to the (n— l)th term 

4tt 2 x 2 ?rV\ 2n) 

very nearly. 

Let m be the greatest integer contained in xw so that d?7r=m-f/> / 
being a proper fraction, then the ratio of the mth to the (m— l)th term 

m 2 /. 3 



* 



(m -\-f)\ 2 m. 

= i -V+ 3 7 

2m 
which is always less than unity. 

The ratio of the (m+l)th to the 'mth term is found in a similar manner 

to be 

2/ 

my 

2 s2 



i+if^/y 
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which is greater or less than unity according as /is less or greater than -^ 
The ratio of the (m+2)th to the (m+ l)th term 

w\4 / 

which is always greater than unity. 

Thus the mth or (m+ l)th term is the least according as / is less or 

greater than -. 
b 4 

When m is large the value of the mth term is very nearly 

1 1 . 2.3.... 2m 

which, by use of the theorem 

1.2.3...*?= ^(2^v)afe" 3e ( 1 + — V 

\ 12#/ 



becomes 



1 . 2s/ir(2m) 2m e' 



Vm ' (2wx) 2m \ 24m J 

=2 \/S-(^)X 1+ 23Si) (i ° 

The number which forms the negative characteristic of the common 
logarithm of this quantity (the mantissa being made positive), reduced by 
unity, denotes the number of ciphers which precede the first significant 
figure in its value ; so that if the mantissa be not made positive, the cha- 
racteristic indicates the number of decimals directly obtainable from the 
series. 

The logarithm of (ii) with its sign changed, after some expansions and 
reductions, becomes 

1 / / / 2 1 \ 
2iU7r#-}-~-lo£; 1n #~- loa-2 — ul -J— + tL_4-_ ] 

f.i being the modulus 43429448 .... Neglecting the last three terms, 
which can only very rarety lead to an error of a unit, we obtain as a result 
that the number of decimal places which the formula (i) will afford directly 
for a given value of x is equal to the greatest integer contained in 

2^7ra?+-log 10 a?— log 10 2 (hi) 

The corresponding value of the least term is, of course, easily obtained 
from (iii), or it follows at once from (ii) ; for the least term 



V w v»+/) 



2,n 
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V(=) 



6 



-2»» 



1 + 



/ 



2m 



m' 



-2h>« 



e -2m~2f 



2 . 

agreeing with (hi). 

The following Table, to replace that in vol, xv. p. 432, was calculated in 
the way indicated above. 



x. 

1 

2 

5 

10 

20 

50 

100 

200 

500 

1000 



Least term. 


Number of decimals 
directly obtainable 
from the formula. 


3rd 


2 


7th 


5 


16th .... 


13 




27 


63rd 


54 


157th 


136 


314th 


273 


629th 


546 


1571st .... 


1365 


3142nd ..... 


2729 



The number of decimal places practically obtainable is limited by the 
difficulty of calculating the Bernoulli's numbers, of which only thirty-one 
have been hitherto obtained. By means of these, however, 156 decimals 
could be obtained of y when #= 100 ; and by deducing a few of the subse- 
quent terms, each from its predecessor (knowing their ratio), 20 places 
more could be obtained without difficulty. 

It is clear therefore that Mr. Shanks' s values of y obtained from #=500 
and #=1000 ought to agree beyond the 59th decimal, if correctly calcu- 
lated. The author at first supposed that the want of agreement was due 
to an insufficient number of the terms involving the Bernoulli's numbers 
having been included, and he therefore undertook the calculation of this 
portion of the expression for #=500 and #=1000 to 100 decimal places ; 
the results, however, still showed a difference in the 59th place. 

In order to determine where the error existed, the same portion of the 

constant was calculated also to 100 places, both from #=100 and #=200 ; 

all the calculations were performed wholly in duplicate, and so much care 

was taken that the author felt a strong conviction of their accuracy, it 

should be noticed that the agreement of all the four results would not 

necessarily prove the accuracy of that value of y ; for any error made in 

11 
the calculation of the harmonic series 1 +- ... . 4- -r^-z would merely pro- 



100 
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duce the same error in 1+- . . . +^77^5 1+tt • • • +7777,, and 1+- . . . 

2 200 2 500 2 

+ -— — , supposing the calculation of the reciprocals beyond -~— to 

have been accurately performed. 

It was therefore evident that the only means of ensuring freedom from 
error in the harmonic series was to recalculate it. The lowest value of x 
which would suffice for a verification was 100 ; the author therefore calcu- 
lated the value of the series 1+- . . . +77^ to 100 places of decimals, 

2 100 

and the result was found to agree to that extent with that given by Mr. 
Shanks in the paper previously referred to ; the value of 1 + 5 ... 4 



2 '50 

was also found to be correct. 

It may be mentioned that the calculation was abbreviated by a simple 
artifice, suggested by Oettinger (Crelle's Journal, t. lx. p. 376), which will 
be easily understood from an example. Suppose the sums of the recipro- 
cals of the odd and even numbers up to 50 are known, and it is required 
to find the sum of the reciprocals up to 100. 

Let 

«=!++-_ . . . 



3 5 49 



then 



2 4 6 50 



^/M+Ui .+!+ l 



so that 



2 2 4 6" 93 100' 

H-i + - . • • + — + — 
2 3 99 100 

_ci4-/3 ,,1,1 , 1 

2 51 53 99 

and the calculation of the reciprocals of the even numbers between 50 and 
100 is rendered unnecessary. 

The values of log 2 and log 5, which were required to form log 100, 
log 200, log 500, and log 1000 were taken from Mr. Shanks' s 'Rectifi- 
cation of the Circle ' (they are also given in the Proc. Roy. Soc. vol. vi. 
p. 397), and the summations of the harmonic series from Mr. Shanks's 
paper in vol. xv. p. 431 of the Proc. Roy. Soc. The results were as 
follows : — - 



#= 100:-— 










y = -57721 


56649 


01532 


86060 


65120 


90082 


40243 


10421 


59335 


93992 


85988 


05770 


42799 


90853 


75858 


22362 


13134 


84454 


84292 


91195 



(A). 
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From x = 


200:— 












r 


= -57721 


56649 


01532 


86060 


65120 






90082 


40243 


10421 


59335 


93995 






35988 


05770 


91390 


77706 


70283 






52607 


02838 


27101 


38332 


64367 . . 


• (B). 


From x = 


500 :— 












y- 


= -57721 


56649 


01532 


86060 


65120 ] 






90082 


40243 


10421 


59335 


93995 






35988 


05771 


53861 


75089 


05973 






20064 


75314 


10504 


07812 


29873 . . 


. (0). 


From x = 


1000:— 












r : 


= -57721 


56649 


01532 


86060 


65120 






90082 


40243 


10421 


59335 


93995 






35988 


05772 


02455 


61942 


00398 






50309 


65017 


53150 


61852 


03044 . . 


• P). 



The terms involving B 39 , B 3l , B 22 , and B ]9 were the highest used in the 
calculation of (A), (B), (C), and (D) respectively. 

It will be seen that (A), (B), (C), and (D) differ in two respects : (i) the 
50th figure in (A) is 2, while in (B), (C), and (D) it is 5, and (ii) all four 
values are totally different after the 59th figure. 

The first discrepancy (i) pointed to an error in the summation of the 

harmonic series. As the author had verified Mr. Shanks' s value of 

1 1 

1 + o • . • + TTwv it was practically certain that (A) was the correct value. 

To place this beyond all doubt, however, the author calculated y from os 
= 50 to 57 places of decimals, and the result entirely confirmed (A). It 
follows, therefore, that Mr. Shanks has made an error of 3 in the 50th 

place in the calculation Taa + TaT • • • + 9oo» so ^at we have the fol- 
lowing errata (vol. xv. p. 431) : — In I + -5- . . . +oa7j> the tenth group of 

five should be . . . 53024 . . . instead of . . . 53027 . . . 

1 1 

In 1 +-7J- . . . + k7Ta the tenth group should be . . . 39677 . . . instead 

of . . . 39680 . . . 

1 1 

In 1 +"9" • • • + iTwIo ^ ie tenth group should be . , . 70880 . . . instead 

of . . . 70883 . . . 

With regard to the second error (ii), which causes the disagreement of 
all the figures after the 59th, it is clear that no error in the Bernoulli's 
numbers could produce such a discrepancy, as the terms involving the 
latter vary their position in each calculation, so that an error in any one 
of them could not affect the same decimal in (A), (B), (C), and (D). 
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By subtraction we find : 

(B) - (A) = . . . (60 ciphers) . 48590 86852 
94425 30244 89703 42646 
54039 73172 (E) 

(D) - (C) = . . . (60 ciphers) . . . 48590 86852 
94425 30244 89703 42646 
54039 73171 . . . (E) 

which only differ by a unit in the 100th decimal place. 

Leaving out of consideration the terms involving Bernoulli's numbers, 
the difference between the series when # = 200 and when #=100 is 

+ T • " + Wo~~ h % 200 )- + 1 ' • ' +ilo - l0 - 10 °) 

11 l 

Similarly, the difference between the series when x = 1000 and when 

x 5= 500 is 

1_ J_ 1 

~~501 + 502 ' ' ' + 1000 g 

Now as it is extremely improbable that two errors, exactly equal in 
amount, should have been made in the calculations of 

11 111 1 

+ -HT3 •••-{- H7T7T and PTTj +F7^ • • • + 



101 ] 102 ' 200 "501 ' 502 n 1000' 

we have very strong evidence that the value of log 2 is inaccurate, and 
that (E) is the correction to be applied to it. 
By subtracting (C) from (A), we obtain 



~(A)=.. 


. (59 ciphers) ... 1 


11064 


84235 


30114 


97702 62179 


26049 


23519 


38678 


*0*»00«0 


• • • • 


• « 



• (F) 

The difference of the corresponding series (omitting as before the terms 
involving the Bernoulli's numbers) 

11 1 . t 



""101 ' 102 1 500 

Having only this one* difference-result involving log 5, it is impossible 
to decide from (A), (B), (C), and (D) whether the harmonic series or log 
5 or both are in error ; but the following reasoning places it beyond all doubt 
that (F) is a correction to log 5, and that the sum of the harmonic series is 
correct. 

1 1 

* By subtracting (B) from (D) we might get another, but the portion oTyr+oAo 

1 

of the harmonic series, as well as log 5, would be common to both,— June 16, 
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Mr. Shanks computed log 2 and log 5 from formulae of the form 

Log 2 = 2 (7 P + 5Q + 3R), 
Log5 = 2(16P + 12Q + 7R), 

P, Q, and R being infinite series (Proc, Roy. Soc. vol. vi. p. 397; ' Rectifica- 
tion of the Circle/ p. 88). 

If, therefore, any error was made in the calculation of P say, it would 
produce errors in log 2 and log 5 proportional to 7 and 16. On trial it 
was found that sixteen times (E) was equal to seven times (F), the difference 
being only such as an error of a unit in the 100th decimal of (E) or (F) would 
produce. This afforded a moral proof that Mr. Shanks had made an 
error equal to one-fourteenth of (E) in the calculation of P, or 



1 



+ 



1 



1 



6 *T" * • • 



31 ' 3.31 3 ' 5.31 

which has rendered his values of log 2 and log 5 incorrect, and that (with 
the exception of the error previously noticed) the harmonic series was 
summed correctly. 

Since log 3 was calculated from the formula 

Log 3 = 2 (11 P + 8 Q + 5 R), 

its value is also incorrect, as also is that of log 10 (log 2 -flog 5) and the 
modulus (the reciprocal of log 10). 

The values of all these quantities are given to 205 decimal places in 
vol. vi. Proceedings of the Royal Society, p. 397 ; but all the figures after 
the 59th decimal place are incorrect in each case. 

The correct values to 100 places are :— 



Log 2 =-69314 


71805 


59945 


30941 


72321 


21458 


17656 


80755 


00134 


36025 


52541 


20680 


00949 


33936 


21969 


69471 


56058 


63326 


99641 


86875 . . . 


Log 3 = 1-09861 


22886 


68109 


69139 


52452 


36922 


52570 


46474 


90557 


82274 


94517 


34694 


33363 


74942 


93218 


60896 


68736 


15754 


81373 


20888 . . . 


Log 5 = 1-60943 


79124 


34100 


37460 


07593 


33226 


18763 


95256 


01354 


26851 


77219 


12647 


89147 


41789 


87707. 


65776 


46301 


33878 


09317 


96108 . . . 


Log 10=2-30258 


50929 


94045 


68401 


79914 


54684 


36420 


76011 


01488 


62877 


29760 


33327 


90096 


75726 


09677 


35248 


02359 


97205 


08959 


82982 . . . 



It is to be observed that the above value of log 3 is not quite as well 
determined as the others ; the calculations in regard to Euler's constant 
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form a real verification of log 2, log 5, and log 10 ; they also verify P, 
Q, and R ; but an error in log 3 in the transcription of P, Q, and R, or 
their multiplication by 11, 8, and 5, or in the final addition and multipli- 
cation by 2, would not be detected. 

In the above logarithms the last figure may be in error to the extent of 
one or two units. 

The value of Eider's constant to 100 decimal places is — 

y= -57721 56649 01532 86060 65120 

90082 40243 10421 59335 93992 

35988 05767 23488 48677 26777 

66467 09369 47063 29174 67497 . . . 

The last figure here also may be in error to the extent of one or two 
units. 

It will be observed that Mr. Shanks' s value* of y for x = 500 differs 

from (C) in the 65th decimal place, and that his value for #=1000 differs 

from (D) in the 73rd place. This is caused by an inaccurate value of B ]3 

8553103 
having been made use of. The correct value of B 13 is g ; but Euler, 

in , , i • -, • 8553103 , A 
who first calculated it, made it ?> " Acta Petropohtana for 1781, 

p. 46), and this incorrect value is given in the ( Penny Cyclopaedia ' 
(Article " Numbers of Bernoulli") and probably in other places. 

The values of the first thirty-one Bernoulli's numbers are given in a 
paper by Ohm (Crelle's Journal, t. xx. p. 11), and B 13 is given correctly 
there. The agreement of the values of Euler's constant contained in this 
paper (when the logarithms of 2 and 5 are corrected) afford a complete 
verification of the Bernoulli's numbers as far as B 26 , and partial verifica- 
tions of the rest. 

The difficulty and inconvenience of making calculations to so many de- 
cimal places is sufficient to warrant the publication of the values of the 
positive and negative parts of the portion of the series involving the Ber- 
noulli^ numbers, in case any one should desire to repeat any part of the 
calculation. We have 

,,■11 1 1 1 , B, B 2 B 3 

' 2 3 x ° 2a? 2x l 4ar ' 6# 6 

and if m denote the sum of the terms of the same sign as the harmonic 
series, and n the sum of the terms of the same sign as the logarithms, 
viz. if 

B x B,, B 5 

2o? 2 bx b 10# ° 

1 B 2 B 4 

2% 4x Scc 8 

* Proc. Eoy. Soo. vol. xv. p. 432. 
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then, when x- 


= 100: — 










m= 


: -00000 


83333 


33337 


30158 


73773 




44885 


67821 


37321 


67823 


08773 




00082 


30639 


33761 


49846 


36254 




14224 


82920 


15089 


52013 


00966 




589 .. . 










n= 


= -00500 


00000 


83333 


33375 


00000 




21092 


80052 


53967 


77636 


59871 




21105 


34750 


25738 


21912 


47654 




32311 


34887 


61384 


97230 


71640 




J-OO « • i 










When* =200: 


'•_—_* 










7)1 — 


: -00000 


20833 


33333 


39533 


73016 




61283 


59538 


59711 


91078 


12258 




22890 


92137 


92329 


17052 


51095 




23490 


89894 


78957 


35141 


06774 




044 .. . 










n~ 


= -00250 


00000 


05208 


33333 


49609 




37505 


14960 


84887 


28877 


09510 




05685 


64006 


82170 


11770 


13023 




01074 


35246 


38458 


75084 


75675 




042 . . . 










When x = 500: 












m = 


: -00000 


03333 


33333 


33358 


73015 




87309 


34487 


73462 


42678 


21147 




87864 


83345 


53789 


61875 


16264 




59799 


20761 


22352 


44796 


75809 




196 .. . 










71 = 


= -00100 


00000 


00133 


33333 


83344 




00000 


00008 


63960 


92825 


14227 




06303 


41363 


50645 


63522 


06215 




12108 


84946 


96404 


61971 


02278 




386 . . . 










When x =1000 


»: — 










m— 


: -00000 


00833 


33333 


33333 


73015 




87301 


59487 


73448 


77428 


21067 




82137 


32165 


00876 


22580 


99464 




52882 


37182 


41954 


63931 


48991 




922 .. . 










n= 


= •00050 


00000 


00008 


33333 


83333 




37500 


00000 


00210 


92796 


09328 




93526 


00040 


82093 


23718 


92820 




00915 


30332 


37713 


49061 


83360 




485 . . . 











The calculation was performed to 105 places, and the last two figures 
have been rejected. 
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Postscript. 

Received June 14, 1871. 

After the completion of the above paper, the author found that Mr. 
Shanks had, in second, third, and fourth supplementary papers on the 
Constant*, extended his calculations so as to determine y from <r=2000, 
5000, and 10,000. 

The values so obtained all differ in the sixtieth decimal ; in fact the 
higher x is taken, the further from the truth are the results, as the errors 
in the logarithms are multiplied by larger factors. 

The calculation for a?=2000 affords a verification of the error of log 2 ; 
for on subtracting the value of y (#=1000) from y (#=2000), we obtain 
(after correcting B 13 ) a result agreeing* with E to 80 decimal places 
(which is as far as Mr. Shanks has calculated the latter value of y), with 
the exception of a difference of a unit in the 73rd figure — an error pro- 
bably in the summation of the harmonic series from 1000 to 2000. 

The values for #= 5000 and # = 10,000 are (besides the errors pre- 
viously noticed) inaccurate from the 62nd figure. 

VIII. " Records of the Magnetic Observations at the Kew Observatory. 
No. IV. — Analysis of the principal Disturbances shown by the 
Horizontal and Vertical Force Magnetometers of the Kew Ob- 
servatory from 1859 to 1864." By General Sir Edward Sa- 
bine, K.C.B., President. Received June 15, 1871. 

(xibstract.) 

This paper exhibits an analysis of the principal disturbances recorded 
by the horizontal and vertical force self-recording magnetometers of the 
Kew Observatory in the years 1859 to 1864, showing the progressive 
diminution in the number and value of the disturbances from a maximum 
in 1859 to a minimum in 1863, being the first moiety of the "decennial 
period ;" and showing also the distribution of the disturbances, increasing or 
diminishing the respective forces, in the several years, months, and hours. 

IX. " Amended Rule for working out Sumner's Method of finding 
a Ship's Place." By Prof. Sir William Thomson, F.R.S. 
Received June 15, 1871. 

In my previous communication on this subject (a?itea, p. 259) I 
described a plan according to which, in the first place, two auxiliary 
lines were to be drawn on the chart, from two sets of numbers taken 
out of a proposed Table, and then Sumner's line (the line on which the 
observation shows the ship to be) was to be interpolated, dividing the 
space between them in the proportion of the differences of the sun's decli- 

* Proc. Hoy. Soc. vol. xvi. pp. 154, 299, vol. xviii. p. 49. 



